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We study the helicity and chirality transitions of a high-energy neutrino propagating in a 
Schwarzschild space-time background. Using both traditional Schwarzschild and isotropic spherical 
co-ordinates, we derive an ultrarelativistic approximation of the Dirac Hamiltonian to first-order in 
the neutrino's rest mass, via a generalization of the Cini-Touschek transformation that incorporates 
non-inertial frame effects due to the noncommutative nature of the momentum states in curvilinear 
co-ordinates. Under general conditions, we show that neutrino's helicity is not a constant of the 
motion in the massless limit due to space-time curvature, while the chirality transition rate still 
retains an overall dependence on mass. We show that the chirality transition rate generally depends 
on the zeroth-order component of the neutrino's helicity transition rate under the Cini-Touschek 
transformation. It is also shown that the chiral current for high-energy neutrinos is altered by cor- 
rections due to curvature and frame-dependent effects, but should have no significant bearing on 
the chiral anomaly in curved space-time. We determine the upper bound for helicity and chirality 
transitions near the event horizon of a black hole. The special case of a weak-field approximation is 
also considered, which includes the gravitational analogue of Berry's phase first proposed by Cai and 
Papini. Finally, we propose a method for estimating the absolute neutrino mass and the number of 
right-handed chiral states from the expectation values of the helicity and chirality transition rates 
in the weak-field limit. 

PACS numbers: 13.15.-|-g, 14.60.Lm, 04.20.-q, 04.70.Bw 



I. INTRODUCTION 

Recent experiments on solar neutrinos [1, 2] have provided strong evidence that flavour osciUations occur while in 
transit to detectors on Earth, and suggest that neutrinos have small but well-defined rest masses. This suggestion 
implies that neutrinos can potentially change their helicity and chirality due to interactions with external fields. It is 
therefore possible to suppose that left-handed neutrinos can be converted to the right-handed form, rendering them 
sterile to the weak interaction, and therefore blind to all known forces except gravitation. 

Although a successful quantum theory of gravity currently does not exist, for energy regimes and space-time length 
scales where quantum mechanics and general relativity can safely appear together, the study of neutrinos in classical 
gravitational backgrounds has the potential to yield, in principle, observationally verifiable predictions. For example, 
their helicity and chirality properties can act as probes of space-time curvature effects appearing as quantum phase 
shifts to identify alternate theories of gravity [3], while other neutrino properties may detect evidence of quantum 
violations of the equivalence principle [4], such as from interactions with matter via the Mikheyev-Smirnov-Wolfenstein 
(MSW) mechanism [5]. In addition, different aspects of neutrino dynamics may reveal valuable information about the 
internal structure of neutron stars and supernovae [6] , and the dark matter content of the known Universe [7] . Because 
neutrinos are produced in abundance by stars and more exotic astrophysical objects like neutron stars-whose surface 
gravity is especially strong-it is important to have a better understanding of their spin and handedness properties 
due to gravitational fields. 

The purpose of this paper is to examine the helicity and chirality dynamics of a high-energy massive neutrino 
while propagating in curved space-time. A previous study [8] of chiral transitions considered the special case of weak 
gravitational fields. It concluded that the transition amplitude is very small in a generally weak background, but 
suggested that it may become large if the background space-time is strongly curved. In this paper, we consider the 
same problem in the presence of strong gravitational fields, but restrict attention to dynamics in a Schwarzschild space- 
time background. An important feature of this approach is to use an ultrarelativistic form of the Dirac Hamiltonian 
derived by the Cini-Touschek unitary transformation [9], which reasonably assumes that the neutrino rest mass is 
very small compared to its total energy. When applied to a free particle in flat space-time, the corresponding Dirac 
Hamiltonian is block diagonal in the chiral representation, where the mass-dependence is contained in the energy 
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eigenvalues, described in special relativistic form. When applied to more general space-times, the mass terms in 
the Hamiltonian can be treated as small first-order perturbations of the massless Dirac Hamiltonian, where we can 
effectively control the degree of mass-dependence as required for a given problem. 

We begin in §2 with a brief review of the formalism behind the Dirac equation in curved space-time that leads to 
the corresponding Dirac Hamiltonian in general form, followed by explicit expressions for Schwarzschild space-time as 
expressed in both regular and isotropic spherical co-ordinates. In §3, we use the derived Hamiltonians to determine the 
helicity and chirality transition rates for an arbitrary fermion, to be later compared with equivalent expressions for the 
high-energy approximation upon developing in §4 the formalism of the Cini-Touschek transformation in curvilinear 
co-ordinates. The helicity and chirality transitions for a high-energy neutrino are then determined explicitly in 
§5. A detailed analysis of the results follows in §6, where special consideration is given to the case of a weak- 
field approximation, which includes the gravitational analogue of Berry's phase first introduced by Cai and Papini 
[10]. Furthermore, we examine the potential role played by the Cini-Touschek approximation on chiral symmetry to 
determine its projected effect on the chiral anomaly in curved space-time. We follow with a conclusion in §7, where 
possible future developments of these results are briefly explored. 

For this paper, geometric units of G = c = 1 are assumed throughout, where we adopt —2 signature for the metric. 
General space-time indices are denoted by Greek characters and range from to 3, while purely spatial indices are 
represented by Latin characters and range from 1 to 3. However, when specifying a particular basis frame, the indices 
are labelled according to specific directions in space-time with respect to co-ordinates which define the chosen frame. 



II. DIRAC HAMILTONIAN IN CURVED SPACE-TIME 



There is a considerable history on the study of spin-1/2 fermions in curved space-time [11], where both gravitational 
interactions and inertial forces due to accelerated motion can be introduced in a self-consistent way [12]. The uniquely 
quantum mechanical nature of intrinsic spin angular momentum makes it a valuable probe of quantum system prop- 
erties in the presence of external classical fields, where predicted physical effects occur for non-inertial fermion motion 
in Minkowski space-time [13] that can be related to gravitational effects via the principle of equivalence. Already, a 
number of well-known experiments [14] suggest this correspondence can be effectively applied here, although it must 
be stressed that critical issues of locality for locally accelerated systems [15] may need to be carefully addressed when 
addressing quantum mechanical issues in curved space-time. Nonetheless, there is no compelling reason to not follow 
this line of reasoning at this time. 

A. Covariant Dirac Equation 

The starting point here is the covariant Dirac equation, 

ni ' 

i^^(x)D^ - -\ ij{x) = 0, (1) 

where m is the neutrino rest mass and Dfj, = -|- iF^ is the covariant derivative operator, with the usual 
covariant derivative on index-labelled tensors and F^ the spinor connection. The gamma matrices {7'' (a;)} satisfy the 
conditions {j^^ (x) , ■y'^ (x)} = 2g'^'^{x) and D^j'^ = 0. Using a set of orthonormal tetrads {e^} and basis one-forms 
{e'^} labelled by caratted indices to define a local Minkowski frame satisfying (e'^, ei>) = df^o , the metric is described 
as g = rjjxoe^ (g) e'^ . The general metric tensor and its Minkowski counterpart are related by vierbein sets {e"^} , 
{e''a} satisfying e" = e"^ e'^ and = e^Q. , such that 

^ ~ ^ /3 ' ^^ct ^ V — ^^v 1 (2) 

9i,u = V&0e"^,eK- (3) 



and 



The spinor connection is 



= -ia"/5(x)r„^^ = -l,.«/5F.^.e%, (4) 

where cr"/^ = |[7",7''] are the Minkowski space-time spin matrices, and T^,^^ are Ricci rotation coefficients derived 
from the Cartan equation of differential forms 

de*^ e" A =0. (5) 
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B. Dirac Hamiltonian in Schwcirzschild Space-Time 



By arranging (1) into the Schrodinger form, it follows that the Dirac Hamiltonian in general space-time co-ordinates 

di ' 



is ih Mr tj){x) = Hil){x), where 



H = ' e% [7^ m + e^p (r?^^ - i a^"") {-ihdj + hTj)] + hTo . 



(6) 



We can now consider a derivation of the Hamiltonian for Schwarzschild space-time. For the metric in Schwarzschild 
co-ordinates, the orthonormal basis frame is 



e* = F{r)dt, e'' = 



1/2 

where F{r) = (1 — 2M/r) is the lapse function. The Dirac Hamiltonian (6) is then 

H = F[m/3+(a-p)]+i^(F-l)aV 



1 



■ dr , e = rd9 , e"^ = r sm9d(j> , 



(7) 



— ih — 
r 

where p = —ih^ in spherical co-ordinates is 

d 



F ( l + r-^lnF^/A a'~ + - cot 61 
\ or J 2 



P 



ifi^, P" 
or 



ih d 

"7 de 



p* 



ih d 



r sin 6 dip 



(8) 



(9) 



and the dot product operation denoted in (8) and elsewhere in this paper is {A - B) = A^B^ + A^B^ + A^B^ . 

For Schwarzschild space-time written in isotropic spherical co-ordinates [16], the one- form frame corresponding to 
the metric is 



where R±{f) = 1 ± M/2 f, and 



^-^''^^dt, = {R+{f)f df , = {R+{f)f fde, 



(i?+(r)) r sin 9 dip , 



(10) 



r(r) = 



9 M 
{F{r)f + F{r) + — 
r 



r(r) = (i?+(r)) r. 
Then the Dirac Hamiltonian in isotropic co-ordinates is 



H 



R- 



TnRll3+{cx-p)- — 



l + r^ln {R_Rl) ) + ^ cot 9 a 



(11) 



(12) 



where p is described by 



d 



9 M' 
iF{r)f + F{r) + - 



d_ 
dr ' 



P 



ih d 
T '39 



P* 



ih d 



r sin 9 dip ' 



(13) 



III. HELICITY AND CHIRALITY TRANSITIONS IN SCHWARZSCHILD SPACE-TIME 
A. Non-Inertial Dipole Operator in Curvilinear Co-ordinates 



Having now the Dirac Hamiltonians (8) and (12) in Schwarzschild and isotropic spherical co-ordinates, respectively, 
we can proceed to determine the helicity and chirality transition rates for an arbitrary uncharged fermion of mass 
m. However, before proceeding it is necessary to first recognize that the momentum operators (9) and (13) have co- 
ordinate dependence, which render them non- commutative. That is, for a momentum operator in general curvilinear 
co-ordinates 



(14) 



where A*(w) arc the scale functions corresponding to the given co-ordinate system, then it is generally true that 
[i/h) = N^^ 7^ 0. Given the three-dimensional Levi-Civita alternating symbol e*-'*^ with indices raised and 

lowered by the Minkowski metric, and where e^^^ = 1, we can introduce [17] a new Hermitian vector operator R, 
which we call a non-inertial dipole operator, for reasons that will be clear later in this paper. The components for R 
in general curvilinear co-ordinates are then 



R" 



2n 

= 5\ 



d 



\nX^{u) - 



d 



1 



d 



\^{u) \du 
1 / d 
A2(u) 



-,lnA^(.) 



1 / d 
A3(u) 



InA^(w) 



InA^(w) 



Ai(u) \dv} 



d 



(15) 



where it is obvious that R vanishes identically for a Cartesian co-ordinate system. For the momentum operators (9) 
and (13), it follows from (15) that 



for Schwarzschild co-ordinates, and 



R"^ = 




cot e ,„ 
= P'', 

r 


(16a) 


R' = 




r 


(16b) 


R^ = 




= —P\ 

r 


(16c) 


R' = 




cote 
= — —p^, 

r 


(17a) 


R' = 




= ip^ 

r 


(17b) 


R^ = 


* r f fli 
^[P ,P ] 


1 a 

r 


(17c) 



for isotropic spherical co-ordinates. It is clear from both representations that R induces some form of apparent 
rotational effect with an associated Lie algebra, due to the mixing of radial and angular momentum states. As well, 
for a fermion propagating solely in the radial direction, R vanishes. 



B. Helicity and Chirality Transition Rates 

Having now determined R explicitly for the Schwarzschild and isotropic spherical co-ordinate systems employed 
here, it is straightforward to compute the helicity transition rate of an arbitrary fermion, as defined hy h = {<j ■ p). 
Adopting the notation 



(18) 



where it is understood that (18) only truly represents the itli component of a three-dimensional cross product when 
applied to Cartesian co-ordinates, it follows from the Heisenberg equations of motion that for an observer very far 
from the gravitational source. 



d 



h = F{-m[— \nF] 13 a'' + -{F- I) {cxxpYp' -F[—\nF^]p'' + ^ cote [I- r—\nF]cx 



h 



d 



dr 



h 



2r^ 



d 



dr 



2i 
r 



d 



dr 



F+{F-1) r— In F^/^ ] {a x pY + - cot 6 {a x p) 



+ i (F - 1) [(a • R) - a'^RH 
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for Schwarzschild co-ordinates, and 
R 



h = 



^^-mi?^ I In (I 



2r 



cot 6* 



2i 
f 



d \ - 1 

1 + f — In R\\{otx pY + - cot 6 {ax p) 



1 - f- In {R.R\) 



1 



(20) 



for isotropic spherical co-ordinates. 

Focussing on the Hcrniitian components of (19) and (20), we notice a paraUel set of contributions to the helicity 
transition rate in botli Schwarzschild and isotropic spherical co-ordinates. There is a noticable exception due to 
presence of the second term in (19) not found in (20). However, this term is non-zero only for propagation that is not 
purely radial with respect to the spatial origin of the space-time background. When this is taken into account, the 
correspondence between the two expressions is more complete. Apart from the co-ordinate singularity for = and 
= TT, the helicity transition rate is finite for fermion propagation in vacuum away from a regular non-rotating star. 
For the special case of a neutron star or black hole, it turns out that while the transition rate becomes singular for a 
fermion near the event horizon in Schwarzschild co-ordinates due to a gravitational blueshift effect, we obtain a finite 
expression in isotropic co-ordinates, where 



iM/2 



M 



1 

4M 



4M2 



cot 6 of . 



(21) 



Given the M^^ dependence on (21), a helicity transition for an arbitrary fermion is essentially nonexistent for 
astrophysically large black holes and neutron stars, but may be significant for microscopic black holes. 

On a more theoretical level, (19) and (20) make an important suggestion that, in the limit as m ^ 0, the fermion's 
helicity transition rate is not a constant of the motion, due to its interaction with the space-time metric during 
propagation. This claim has also been put forward [17] in the context of accelerated fermion beams in circular storage 
rings, where the external magnetic field required to bend the beam generates the violation of helicity conservation. 
Therefore, it should not be too surprising to obtain a similar result due to a non-trivial classical space-time background. 
In sharp contrast, the chirality transition rate is a constant of the motion in the massless limit, since for the projection 
operator P± = (l ± 7^) /2 where the upper sign refers to the right-handed state. 



(Schwarzschild) 

(Isotropic) 



which has an overall mass dependence. Furthermore, it is known that for [8] 

\i>c) = Cr|(^r) + Cl|<^l) , 

where 











(22a) 

(22b) 

(23) 

(24) 



are the right- and left-handed chiral states subject to the normalization condition |Cr|^ -|- |Cl|^ = 1, it follows that 
(/37^) = (i/'cl Pl^ IV'c) = 0. That is, the chirality remains a constant of the motion, even for massive fermions. This 
discrepancy challenges the well-known identification [18] between helicity and chirality for massless fermions, and 
merits further investigation. 



IV. CINI-TOUSCHEK TRANSFORMATION AND THE HIGH-ENERGY DIRAC HAMILTONIAN IN 

SCHWARZSCHILD SPACE-TIME 

Although (19) and (20) are obviously useful as exact expressions for the helicity transition rate of an arbitrary 
uncharged fermion, it is also useful to reconsider the problem in terms of a high-energy approximation of the Dirac 
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Hamiltonian. which can be written as a power series expansion with respect to q = m/p <^ I. By doing this, we can 
identify the leading-order terms that contribute to the dynamics of particles which satisfy this condition. High-energy 
massive neutrinos certainly qualify for this treatment. In addition, by adopting the chiral representation [18] 




we seek to re-express the Hamiltonian such that the zeroth-order term is in block diagonal form (while still retaining a 
particle mass dependence), where the first-order term is a mass perturbation. This allows for an effective decoupling 
of the four-spinor wavefunction into its right- and left-handed states, while the perturbation term is responsible for 
their mixing. It is known that, by applying the Cini-Touschek (CT) transformation to the free-particle Hamiltonian 
[9], the new Hamiltonian effectively removes the /3-matrix term and yields an energy eigenvalue E = ±y/p'^ + 'w?. 

In this Section, we re-derive the CT transformation for the free-particle Hamiltonian in curvilinear co-ordinates, 
then use it to evaluate the high-energy form of the Hamiltonian in Schwarzschild space-time and its weak-field 
limit. Unlike the more well-known Foldy-Wouthuysen (FW) transformation [19], used extensively to model the 
non-relativistic motion of a spin-1/2 particle in external fields, the CT transformation has not been as widely used 
for similar purposes. A preliminary study has recently been considered [20] in the context of non-inertial motion 
through Minkowski space-time, where the CT transformation is performed on a local inertial tangent space in a 
frame comoving with the particle, with the outcome described by general co-ordinates after projecting back onto 
the space-time manifold. The approach taken here instead follows the known procedure when applied to the FW 
transformation, as adopted earlier [8] in the study of neutrinos in weak gravitational fields. 



A. Free-Particle High-Energy Hamiltonian in Curvilinear Co-ordinates 



To derive the CT transformation for the free-particle Dirac Hamiltonian [17], we begin with the unitary operator 
exp {iSct), where 



i>CT = —z — P{a-P) 
2p 

and u){q) is a constraint function used to absorb the /? matrix. Then 



HcT = e'^c-^ [m/3 + [ol ■ p)] e 



[m/3 + (a • v)\ ■ 



By Taylor expansion. 



-/3(a-p) 

V 



1 
3! 



-/3(a-p) 
p 



V 



■fi{a-p) 



[{p-p) + h{a-R)] . 



(26) 



(27) 



(28) 



(29) 



With this expression, cr ■ R la analogous to a magnetic dipole term due to the curl of the electromagnetic vector 
potential, hence the reference to i? as a non-inertial dipole operator. 
By substituting (29) into (28), we show that 



cos 



w 1 + 



2p2 



{a-R) 



sm 



2p2 



{CT-R) 



1 - 



{OL-p) 



Applying (30) to (27) and setting the coefficient of /3 to zero, it follows that 

h 



>{d) « ( 1 - ^ («■ • ^) ) tan-i 



2p2 



with the final result that 

e'-^cT [m(3 + [a ■ p)] e''^^-^ 



\/TTg2 2m 



(a-R) 



(a-p) 



(30) 



(31) 



(32) 
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For the small angle approximation of (31), 

uj{q) ^ q(^l-^ {a-R)y (33) 

In Minkowski space-time, the presence of the dipole term in (32) suggests that a small energy shift appears for 
a spin-1/2 particle in the ultrarelativistic approximation. For high-energy neutrinos, this shift should be negligible 

compared to its rclativistic energy ^/p^ + m?, since the dipole term is order q'^ <C<C 1. However, for neutrinos with 
lower energies and other spin-1/2 particles with larger rest masses, it may be a detectable effect for large enough q. 



B. High-Energy Hamiltonian 

To derive the high-energy form of some operator using the CT transformation, we use the series expansion [19] 



*2 '3 

XcT = e'^^-Xe-'^^- = X + i[ScT,X] + '-[ScTAScT,X]] + '-[ScT,[ScT,[ScT,X]]] + ' 



(34) 



where Sct is of order q and X is an arbitrary operator. Because the CT transformation operator is unitary, the 
eigenvalues associated with the transformed operator remain unchanged [21]. For q sufficiently small, as satisfied by 
the case for massive neutrinos, it is sufficient to include only the first-order expansion term in (34). Therefore, 



w(g) 
2p 



[/3{oc-p),X]. 



(35) 



We therefore insert 1 = e^^'^'^ g-^-ScT between each factor and operator that comprises the given Hamiltonian and 
evaluate to first-order in uj{q) using (35), keeping terms only up to first-order in H. 
By applying (32), (33), and (35) for X replaced by H, it follows that 



HcT 



) + (i^-l)a>''-- 



F ( l + r-^\iiF^/A a'' + 1 cotOa 



P 



-(F-l)(l 

9" 



{a-R)] p^p 



dr 

h 

2 



dr 



InF (o- X pY -{F-l) [{a ■ R) - a'^R^ 



+ ^ 9^ (a • P) - „ 

2 p 2m 



for Schwarzschild co-ordinates, and 



R- { , . ih 



[7' {R-p)+icx-{Rxp)]^ 



(36) 



1+f In (R_RIY'" a^ + -cotea^ 
or ^ J 2 



Ri 



{Rl-l)ll--{a-R)]{p-p) + h 



d (R^ 



1/2 



(cr X pY + {a-R) 



-q^ {a-p)-^—[-f^R-p)+ia-{Rxp) 
2 p 2m 



(37) 



for isotropic spherical co-ordinates. As expected, we retain the zeroth-order part of the original Hamiltonian, while 
generating the first-order mass perturbation terms proportional to (3, which couple the right- and left-handed states of 
the four-spinor. In addition, we have higher-order mass corrections to the diagonal parts of the Hamiltonian, including 
a contribution due to R. For neutrino interactions, we neglect the higher-order mass terms in (36) and (37). 

V. HIGH-ENERGY HELICITY AND CHIRALITY TRANSITIONS IN SCHWARZSCHILD 

SPACE-TIME 



In order to determine the high-energy form of the helicity and chirality transition rates, we need to first compute 
the helicity and projection operators in the high-energy approximation from (35), since 



^CT = ^[Hct,Xct] = e'^°--[if,X]e-'^^ 



(38) 
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Then it follows that to first-order in ij{q), 

hcT « h, 



2 p 2 p 

r(0) , oj(q) 



If we write the high-energy Hamiltonians (36) and (37) in the form Hct = H^t + ^ /^^ct ' tlien 

m— P fl 

h 



Pc 



(39) 
(40) 



(41) 
(42) 



It is interesting to note from (40) and (42) that the chirality transition rate in the high-energy approximation is 
dependent on the zeroth-order contribution of the helicity transition rate, the implications of which are unclear at 
present. 

By straightforward application of (38) with (39), we show that the leading order expressions for the high-energy 
helicity transition rate for a neutrino in Schwarzschild space-time, as seen by a distant observer, are 



TO=0 p 



hcT « h 
in Schwarzschild co-ordinates, and 
hcT ~ /I +- -53 

m=0 p 



+ 1 (F-^\nF]{l-i^{a-R)]f3[{2F-l){a^p^)p^ + {a-p)p^-a^p-p) 



dr 



■In 



R^ 



{a-p)p^-Rlpn 



R- 

r; 



o-'' (p ■ p) 



(43) 



(44) 



in isotropic spherical co-ordinates. In both (43) and (44), we notice that the high-energy corrections to the zeroth- 
order helicity transition rate arc dependent on the helicity operator and projected in the radial direction. As well, 
these corrections reduce to zero in the limit as M ^ 0, as we should expect. For a neutrino near the black hole's event 
horizon, we find that the high-energy correction for (43) also grows singular due to gravitational blueshift, while for 
(44) it is shown that 



hcT 



fRiM/2 



1 

'4M 



l--/3(T-p) 
P 



(45) 



Similar application of (38) with (40) shows that the leading-order Hermitian contributions to the high-energy 
chirality transition rate in the high-energy approximation are 



PcT± « ± - F/3 
P 



l{F-l){axpY p^+rf^ \nF-']p^ + j 



1 + r— In 
or 



9 ^ ^l/2\ 



P 



1 

27 



cot6'p'' 



and 



■PcT± 



P RV 



dr 



In 



R- 



+ f-^ In {R-RlY^') p' + ^_cote p' 
r \ or ^ J 2r 



(46) 



(47) 



in Schwarzschild and isotropic co-ordinates, respectively. Proceeding similarly as with the helicity transition rate, it 
is straightforward to show for both (46) and (47) that 



CT± 



M=0 p 



-p^ + — cot0p' 
r 2r 



in the limit M — > 0. As for neutrino propagation near the black hole's event horizon, it is shown that 

PcT± « ±-i-^/3(7' + 4)p'- 

rRi2M 8M p ^ ' 



(48) 



(49) 
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for Schwarzschild co-ordinates, while 

Pgt± 



« ± - /? (7^ + 2) (50) 



for isotropic spherical co-ordinates. Apart from a factor of two in the second terms of (49) and (50), the chirality 
transition rates at the event horizon are identical. 

There is one difficulty, however, worth mentioning as it concerns isotropic co-ordinates near the event horizon, 
in that the radial momentum operator becomes degenerate when expressed in terms of the Schwarzschild radial co- 
ordinate, as shown in (13). This makes the interpretation of (45) and (50) more difficult because the co-ordinate 
system is ill-defined at the event horizon, once the f-dependence is rewritten in terms of r. 

VI. SPECIAL CONSIDERATIONS 
A. Weak-Field Limit with Gravitational Berry's Phase 

In order to identify neutrino helicity and chirality interactions in an observationally realistic context, it is necessary 
to consider the weak-field limit of curved space-time. To do this requires that we perform a first-order expansion 
about Minkowski space-time, such that the metric is described by g = {rjap + hap) dx" dx^ , and the covariant Dirac 
equation (1) takes an equivalent form for a local tangent space centred about the gravitational source. An interesting 
development proposed by Cai and Papini [10] is that the phase of the wavefunction can contain gravitational field 
information that leads to a covariant generalization of Berry's phase [22] when considering closed particle trajectories, 
with an explicit expression for the case of the weak-field limit. It is further shown [20] that the gravitational phase 
can be introduced into the Dirac Hamiltonian as a vector gauge field within the framework of Minkowski space-time. 

While a more complete treatment of the gravitational Berry's phase is found in Appendix A, it suffices to state 
that (1) is equivalent to 



nMVA + *r^ + ^(VA$G))-f 



V(a;) « 0, (51) 



where $Gj the gravitational Berry's phase, takes on a non-zero value only for closed trajectories in space-time [10]. 
Because $g was originally derived in terms of a Cartesian co-ordinate frame, it is best to consider the weak-field limit 
of this problem in the same fashion. Therefore, given the one-form basis set corresponding to the weak-field limit in 
Cartesian co-ordinates, 

-0 = ri-?M)"^,. 



2M\ 



1/2 

da;*~ , (52) 



where {x^,x'^,x^) = {x,y,z) and r = \J x^ -\- xp' -\- z^, a straightforward calculation shows that the Dirac Hamiltonian 
up to first-order in M/r is 



H « m(l-^^/3+(^l-^^ (a-p) + ^^(a-Vr) + (a- Vc1>g) + (Vo3>g) , (53) 



where = —ihd/dx^. While it is clear from Appendix A that the gravitational Berry's phase is Lorcntz invariant, it 
is important to emphasize here that the gradients of $g sje frame- dependent and will vary in appearance for different 
co-ordinate systems. Applying the CT transformation to (53), we find that 

HcT ^ {l- (a • p) -h y ^ (a • Vr) -h (a • V$g) + (Vq^g) 



P 



[p ■ p) - h-;^ cr ■ ( Vr X p) 



r 



\q-{l-M.ya.-p), (54) 



where the second term in the linear mass-perturbation contribution to (54) indicates a spin-orbital coupling interaction 
present. It follows that the helicity transition rate for a high-energy neutrino in the weak-field limit of Schwarzschild 
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space-time is 



hcT ~ h 



^=0 p 



( Vr ■p){(T-p) (cr • Vr) [p ■ p) 

r r 



(55) 



where 



M 2M 2 r 9 

h = -m -^(3 {a- Vr) - {Vr -p)- -a- ( V$g xp)-T v *G - <t • V (Vq^g) 



o • A'/ ^ 5 ^ 

— 3 1 — a • ( V r X p) + i/i7 — 



- (Vr ■ Vr) - - V^r 



(56) 



As well, the weak-field approximation of the chirality transition rate for a high-energy neutrino is 



PCT± « ±- /3 

p 



2M 

T2~ 



1-77) • P) + ^ " • C^^G X p) + 7^ V-'^G + o- • V (Vo^'g) 



(57) 



For comparison, wc can derive an expression equivalent to (53) in spherical co-ordinates by performing a first-order 
expansion of M/r in (12). Then the Dirac Hamiltonian becomes 



H « m 1- 



M 



2M\ 



r I 



(a • p) + 



1 H OL — 



2r ) 



ih 

2r 



1- 



2M 



cot e a" 



+ (a • V$g) + (Vo$g) 



(58) 



in the weak-field limit, where p is described by (9). It is trivial to note-but important to emphasize-that the free- 
particle Hamiltonian expressions from (53) and (58) are distinct in the limit as M/r 0, due to the choice of 
co-ordinate system. Again, by applying the CT transformation, we obtain for the high-energy approximation of the 
Hamiltonian 



-ffcT « 1 (a-p) + - 1-— a'^-— 1 cot0a^ + (a- V$g) + (Vo$g) 



2r ) 



2r 



r / 



V 



— [{p ■ p) + fi{cr ■ R)] + h— {a X pY 

r r 



M 



o9 



(59) 



where a similar type of spin precession term appears in the linear mass-perturbation correction. Then the high-energy 
helicity transition rate follows as 



hcT ~ h 



1=0 p 



2M M 

{cr p)p cr {p-p) 

r r 



(60) 



where 



h = —m — pcr 



2M ^ h 4M\ 
— P +^'1 I 



r / 



cot ^ a"^ - - a • ( V$G X p) - Y (V • V$g) - cr V (Vq^'g) 



I 
r 

ih r 

-^1 



2 ( 1 - j (a X p)"^ -I- ( 1 — 1 cot 6» (a X p) 



i (a X V) • V$G 



1 



5M\ 1 1 



- + - 



r ) 2 sin^ e 



1 



(61) 



and V ■ V = Vj + V 2V 2 + VgVg , which is only equal to the Laplacian operator in Cartesian co-ordinates. 
Similarly, the expression for the chirality transition rate of a high-energy neutrino is 



CT± 



2M 



P -7 



r \ 2r ) 2r \ r ) 



+ _ a . (V$G X p) 7^ (V • V$g) + o- • V (Vo$g) 



(62) 



A comparison between the two sets of calculations in the weak-field limit suggests that, in the limit as Mjr ^ 0, the 
helicity defined in Cartesian co-ordinates is a constant of the motion, while the same is apparently not true in spherical 
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co-ordinates. However, this is merely a co-ordinate effect that comes from the spinor connection in Minkowski space- 
time. To demonstrate this more precisely, we can calculate the expectation value of the helicity transition rate for a 
local observer comoving with the neutrino, such that 



^ = [500(3:)] h 



(63) 



for proper time r. We can define a normalized state vector \'i'±) = |±) (Xi \ipc) , where lipc) is described by (23), and 



1+) 



-) 



(64) 



are the neutrino's positive and negative helicity states with respect to its quantization axis in the radial direction, 
such that 

(=F|c|±) = [cos COS =F * sin (/?] ^ -|- [cos sin ± i cos (p] y — sin z (65) 
for Cartesian co-ordinates, and 

(t|o-|±) = x''±ix'^ (66) 

for spherical co-ordinates. Then it can be shown that the high-energy helicity transition rate is {dh/dT)Qr^^ = 
(^'q:! {dh/dT)Q,r^ l^±) J where for Cartesian co-ordinates it follows from (55) and (56) that 



/dh\ 

\dr/cT± 



-(T|cr|±) 



+ 2Re{CnCl) 



V(Vo$g; 
2M 



(|Cr|^-|Cl|^) 



^ ( V$G xp) + 3i^ ( Vr X p) 



P 



M 



( Vr ■ p)p Vr (p ■ p) 

r 



(67) 



It is obvious here that (67) vanishes when M/r 
frame, we can show from (60) and (61) that 



0. Performing the same analysis for the spherical co-ordinate 



dh\ 

'^^ / CT± 



-{t\t\±) 

(|Cr|' - |Cl 

2Re {Cb CD 



V(Vo$g) + (|Cr|'-|Cl|') 



(V$G X p) + h{R^G) 



1 



1 - 



M 



cot 6p^ — 



ih 



2r2 



^_3MX 



r / 



cot 6* 



{p'±ip'^)p\ 



(68) 



'±ip'^) - - 
r \ 'AT J " r 

q 2M 

p r 

which is generally non-zero. However, without any loss of generality it is evident that, for exclusively radial propagation 
{p^ = p'^ = 0) along the 6 = 7r/2 plane, it follows that (68) vanishes as well in the M/r — > limit. Similarly, it is 
straightforward to show that the helicity-flip induced chirality transition rates for (57) and (62) are 



!T± 



T2Re {CnCD^ (A ± iy) • V (Vo$g) 

P 



in Cartesian co-ordinates, and 



CT± 



T2Re {CnCl 



)-( 
P ^ 



V(Vo«'g) 



(69) 



(70) 



in spherical co-ordinates. 

A comparison of (67)-(70) shows that the gravitational Berry's phase makes a meaningful contribution to the high- 
energy helicity and chirality transition rates. In particular, the term V (Vo$g) persists regardless of the neutrino's 
degree of handedness or mass dependence where it concerns the helicity transition rate. One intriguing possibility that 
arises here is the potential ability to determine the absolute mass of a neutrino and the amount of right-handedness 
within a typical neutrino beam from measuring the helicity and chirality transition rates, and performing a parameter 
fit on the data to determine Cr , Cl , and q . In practice, however, such experiments would be extremely difficult to 
perform because of the need for an enormous amount of neutrino flux to get statistically significant measurements. 
According to the Standard Model [23], neutrinos have only negative helicity and are exclusively left-handed. Given 
the Superkamiokande and SNO evidence in favour of neutrinos with rest masses, we should expect some right-handed 
states to exist (at least for Dirac neutrinos), but to the degree that |Cr| ^ 1 and |Cl| ^ 1. Unless this assumption is 
shown to be completely wrong, there will be practically no linear mass contribution to the neutrino helicity transition 
rate and effectively no observable chirality transition rate whatsoever. Therefore, it would be a great surprise if we 
were to obtain an unambiguously non-zero measurement of the helicity and chirality transition rates with discernable 
mass dependence, based on (67)-(70). 
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B. Chiral Anomaly 

Although this paper is primarily interested in the quantum mechanics of neutrinos in curved space-time, there are 
legitimate reasons to briefly consider potential implications to quantum field theory in this context. It is beyond the 
scope of this paper to consider a detailed treatment of neutrino quantum field theory in a gravitational background, 
but one particularly relevant issue to examine is the possible impact of the CT transformation on the chiral anomaly 
[24] in curved space-time. It is understood [25] that the conservation of the chiral current in the massless limit is 
violated by the Pontrjagin topological invariant, which is proportional to the square of the Riemann curvature tensor. 
Here, we examine whether the CT transformation has any significant impact on the known expression and discuss its 
possible implications. 

To begin, we have the action for the Dirac field 



where the Lagrangian density £d is 



m 



and tlj{x) = tp^x) (3. We know that if we perform an infinitesimal chiral transformation [26] 

ip{x) cxp [i€{x) 7^] tp{x) PS [l + ie{x) 7^] i^{x) , 
y^^'tpj^ 'y'-^{x)tp satisfies the equation 



the chiral current 







If we now perform the CT transformation on (74), we obtain to leading order in q 
f 



1 „ 2m - 5 

7= 9i,J^CT - V'CT 7 V'CT 



1 



P V-9 



■ da 



-9 tpCT 7 



-/3[(a-V)7l 



e^j (<T X pY -i e^Q (a • p) 



where 



(71) 



(72) 



(73) 



(74) 



(75) 



•^5CT 



i fAcT 7^ 7^ (x) -ipci 



(76) 



We see from (75) that the mass-dependent correction to the chiral current is frame-dependent, due to the presence of 
p in its definition. Clearly, when taking the massless limit, we obtain the conservation of chiral current which follows 
from (74). Therefore, it follows that while there is a small mass correction to the original chiral anomaly equation 
[25], we expect no impact on the chiral anomaly in the massless limit, since the Pontrjagin invariant comes from the 
spinor connection F^, which is not present in (75). 

It should be noted that there is an ambiguity in determining the modified chiral current equation, because if we 
instead performed the CT transformation first, followed by the chiral transformation, we obtain in place of (75) 







; ^i^JsCT - 7^ V'CT 



_ q 

p 

- V'CT 7^ 



1 



-5 V'CT 7 



h 



l3[{oc-V)r]-ie"jP^ 



9a V'CT - it/jCT 7^ 



V'CT 

^/3[(a-V)7"r„] + |{/3a,7«ra}.p 



V'CT 



(77) 



which has a non-trivial dependence on . However, it seems more correct to follow the approach leading to (75) 

because the original chiral current equation (74) follows from an exact Dirac Lagrangian density, while the chiral 
transformation leading to (77) comes from an approximation of £d, so the correspondence is not exact. 
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VII. CONCLUSION 



In this paper, we determined the hehcity and chirahty transition rates for a high-energy massive neutrino propa- 
gating in a Schwarzschild space-time backgronnd, where the dynamics is governed by a Hamihonian transformed by 
the Cini-Touschek unitary operator. We derived the high-energy Hamiltonian to first order in g = m/p in terms of 
both regular Schwarzschild co-ordinates and isotropic spherical co-ordinates, and determined the helicity and chirality 
transition rates accordingly, with particular attention paid to physical effects near the event horizon of a black hole. 
Under general conditions, we learned that the neutrino's helicity transition rate in a gravitational background is not 
a constant of the motion for massless particles, while the chirality transition rate still retains an overall dependence 
on neutrino mass, and thus vanishes in the massless limit. It was demonstrated that a relationship exists between 
the chiral transition rate and the zeroth-order contribution of the helicity transition rate in the high-energy approx- 
imation. We further considered the special case of a weak gravitational field with the inclusion of the gravitational 
analogue of Berry's phase, where the expectation values for the helicity and chirality transition rates were obtained 
and compared. In addition, we determined the mass-dependent corrections to the chiral current equation and its 
relationship to the chiral anomaly when making use of the Cini-Touschek transformation. 

For this investigation, we used the helicity operator as defined in flat space-time. Strictly speaking, the more 
appropriate operator to use is the time component of the Pauli-Lubanski four-vector [18] defined in curved space- 
time. The next development for this research is to compute the covariant spin dynamics of neutrinos in Schwarzschild 
and Kerr backgrounds using the Pauli-Lubanski vector [27], which can be compared to the results contained in this 
paper. Cosmological considerations during the earliest stages of the known Universe may be relevant for applications 
of this type of research. Another possibility of future development is to consider a detailed study of neutrino wave 
packets within this context, to see how the spin and chiral behaviour would diflfer for a wave packet description, as 
opposed to plane waves. These and other types of possibilities may be developed in the near future. 
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APPENDIX A: GRAVITATIONAL BERRY'S PHASE IN THE WEAK-FIELD LIMIT 



Berry's phase [22] was originally formulated in terms of the non-relativistic Schrodingcr equation, assuming adiabatic 
transport of the quantum system around a closed path in parameter space. It is a geometric and gauge-invariant 
phase that is identified with the motion of the system in the space of parameters. One generalization of this formalism 
[28] is to evolve the system using Fermi- Walker transport in place of parallel transport, which has the advantage of 
applicability for non-geodesic motion. 

The covariant generalization of Berry's phase, as developed by Cai and Papini [10], allows for the introduction of 
Lorentz-invariant classical fields within a quantum framework. One example of this is the covariant generalization of 
the Aharonov-Bohm effect [10, 29]. In this Appendix, we present a brief description of the gravitational Berry's phase 
for the weak- field limit of Schwarzschild space-time. This is formally described by the line integral 



^ Jxo ^ Jxo '- 



(A.1) 



where the integration path is specified along the geodesic that describes the particle's motion, p** is the momentum of 
a free particle described by the wavefunction and L"^ is its corresponding orbital angular momentum satisfying 



{x^ - z")p^ - {x^ - z^)p^ 
ihd" *o- 



Formally, it happens that when considering transport of the quantum system in a closed path C in space-time 
according to (A.l), it can be shown thatthe gravitational Berry's phase can be identified with 



(A.2) 
(A.3) 



*g(C) 



a/3 / 



(A.4) 
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where 



a0 



dz^h 



dz'^r 



&0\ 



/3,A + ^SAJ ^/3A,a 



are the translational and rotational components of the Cartan connection [30]. 
In Cartesian co-ordinates, it follows that the most general expression for (A.l) is 



where 



(V^$g) 



M , „ 2M' 



M 



By a simple co-ordinate transformation, the gravitational Berry's phase in spherical co-ordinates is 

px px px px 

$G = / dt' (Vt$G) + / dr' (V^$g) + / de' (V^^g) + / V (V^*g) , 

Jxf\ J Xn J Xn J Xn 



where 



(V«$g) 



--^P* + ^ [cose* cos 9' + sin 6* sm.e' cos(<^ - <^')] [^P* - - ^Op*^] > 



r' 
M 



(V,$g) 

(Ve$G) 

(Vv,$g) 

for p* = —ihd/dt and p'' = —ihd/dr. 



— [cos 9 cos 9' + sin 6 sin ^' cos(y — (p')] . 

2M [cos 9 sin 9' - sin 6 cos 6*' cos((/j - tp')] p"" , 
-2M sin 6* sine*' sin(<^ - ip')p'^ , 



(A.5a) 
(A.5b) 



$G = r dt' (Vt$G) + dx' (V,$g) + rf?/' (Vj,$g) + (V,$g) , (A.6) 

J Xn J Xn J Xn J Xn 



(A.7) 
(A.8) 

(A.9) 

(A.10) 

(A.ll) 
(A.12) 
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